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ABSTRACT 

In a matched sample of Seyfert and quiescent galaxies we simultaneously probe the stellar and 
molecular gas kinematics from 1 kpc down to 50 pc with the aim of identifying the dynamical processes 
dictating black hole accretion rates. This first paper compares the global characteristics of a sample 
of ten galaxies. We find several differences within a radius of 500 pc that are correlated with Actkve 
Galactic Nucleus (AGN) activity. The Seyferts have: (1) a more centrally concentrated nuclear stellar 
surface brightness with lower stellar luminosities beyond a radius of 100 pc, (2) a lower stellar velocity 
dispersion within a radius of 200 pc, (3) elevated H2 1-0 S(l) luminosity out to a radius of at least 250 
pc, and (4) more centrally concentrated H2 surface brightness profiles. These observed differences can 
be interpreted as evidence for Seyfert galaxies having a dynamically cold (in comparison to the bulge) 
nuclear structure composed of a significant gas reservoir and a relatively young stellar population. 
This structure is undetected (and possibly does not exist) in quiescent galaxies. The presence of 
such a nuclear structure in Seyfert galaxies provides evidence for inflow of the surrounding interstellar 
medium since the nuclear stellar population requires a supply of gas from which to form. The fueling 
of a Seyfert AGN is thus associated with the formation of a dynamically cold component of gas and 
stars on scales of hundreds of parsecs. 

Subject headings: galaxies: active — galaxies: kinematics and dynamics — galaxies: nuclei — galaxies: 
Seyfert — infrared: galaxies 



1. INTRODUCTION 

\ We have now come to the understanding that super- 
massive black holes (BHs) play a pivotal role in galaxy 
evolution. This revelation has been driven by the dis- 
covery that most, if not all, galaxies harbor a BH, and 
that the masses of these BHs are r elated to the global 
■ properties of their host galaxy (e.g. iFerrarese h Merritti 
; [200alGebhardt et al.|[2000llTremaine et al.ll2002[ ). This 
relationship is thought to be a consequence of the coevo- 
lution of the BH and its host galaxy, and identifying the 
origin of this coevolution is an essential component to our 
understanding of galaxy formation and evolution. It is 
widely accepted that BHs accumulate their mass through 
an Active Galactic Nuclei (AGN) ph ase in which they 
\ accret e material from the host galaxy. iKaufFmann et al] 
. (|2003| ) find that Seyfert AGN generally reside in host 
galaxies with a younger stellar population relative to that 
found in quiescent galaxies, confirming that an abundant 
fuel supply is available in these galaxies on scales of sev- 
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eral kiloparsecs. However, we have yet to identify the 
processes responsible for actually driving material down 
to the radii necessary for accretion onto the BH. We must 
therefore first identify and understand these processes 
before forming a complete picture of the role of BHs in 
the evolution of galaxies. 

Major mergers have been shown to be an efficient 
mechanism for rapid loss of angular moment that can 
suppo rt the fueling of the most l uminous quasa r AGN 
(e.g. iBarnes fc HernquistI [TOOl iHopkins et al.l [2008) . 
However, it is likely that this merger driven inflow is 
stalled at scales of 1 kpc, at which point it is actu- 
ally smaller scale processes that are expected to be 
responsibl e for dire c tly fu eling the central BH (e.g. 
[Hopkins &: Quataertl l2010f ). In addition, multiwave- 
length surveys find that it is only the most luminous 
AGN that are associated with major mergers leaving a 
significant fraction of AGN out to a redshift of at least z 
= 3 to be driven by other processes (e.g. Treister et al. 
2012). Recent studies have also shown that at redshifts of 
z > 1 a significant fraction of AGN reside in disk galaxies, 
which strongly suggests that secular disk processes rather 
than major mergers are the principle driver fueling their 
AGN (e.g. iSchawinski et al.ll201lL l2012llKocevski et al.l 
120121 ). By redshifts of z ~ 1-1.5 disk galaxies begin to 
show the characteristics of the local Hubble sequence 
(Ivan den Bergh et alj 120001 iKaiisawa fc Yamadal [200l 
iConsehce et al.l 12004 lOesch et al.l I2010D and bars can 
be traced out to z ~ 1 ( Sheth et al.ll2008[ ). thus secu- 
lar processes become a promising alternative source of 
AGN fueling by these redshifts. Although secular pro- 
cesses capable of driving the needed quantities of fuel 
inward have been identified, e.g. large-scale spiral struc- 
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tures and bars, these processes are only thought to op- 
erate down to scales of hundreds of parsecs. This may 
explain why previous studies of large-scale structures in 
local galaxies have f ound little, or no, cor r elation with 
AGN activity (e.g. [Sh losman ct al.' '2000', 'Lain e et all 
120021 [Laurikainen et al.. .2002. 2004, Hao et al. 200"9[). It 
is thus a series of processes that c arry material inward 
(jShlosman et al.lll990llWadall2004D and it is the smaller 
scale processes (less than 1 kpc) that are ultimately reg- 
ulating BH growth in the local universe as well as at high 
redshifts. 

To constrain what processes are relevant to fueling lo- 
cal AGN, and are likely candidates for fueling at higher 
redshifts as well, we have observed a sample of matched 
Seyfert-quiescent galaxy pairs from which we can identify 
properties unique to Seyferts that are potentially related 
to fueling of their AGN activity. A significant impedi- 
ment to addressing this issue thus far has been the limita- 
tion of relying on morphological signatures of these pro- 
cesses, while the critical difference between Seyfert and 
inactive galaxies is expected to lie in the kinematics of 
the interstellar medium (ISM) . Comparison studies of the 
circumnuclear dust structures in the central kpc of galax- 
ies imaged with the Hubble Space Telescope (HST) have 
found little or no correlation between isophotal twists, 
nuclear bars, and nuclear spiral s detected in Seyferts 
relative to inactive gala xies (e.g . 'Hunt &: Malkanl j2004 
IMartini et all [2003a! bV ISim5es L opes ct al. (200'q) find 
no correlation between the presence of circumnuclear 
dust and AGN in late-type galaxies (the typical hosts 
of Seyfert AGN) , but do find a correlation in early-type 
galaxies. However, they conclude that even in early- type 
galaxies this structure is not a sufficient condition for 
AGN activity since similar structure is also seen in a sig- 
nificant fraction of comparable quiescent galaxies, and 
that there must be other critical factors in driving AGN. 

Exte nding beyond morph ological signatures of AGN 
fueling, iDumas et al.l ()2007D carried out a 2D kinematic 
comparative study of a sample of Seyfert and inactive 
galaxies and report that on scales of 200 pc the kinemat- 
ics of the rotating disk of ionized gas in Seyfert galaxies 
deviate from axisymmetric rotation to a greater degree 
than in inactive galaxies. They conclude that the link 
between the AGN and the host galaxy dynamics there- 
fore exists on scales of hundreds to tens of parsecs. This 
suggests that these are the dynamical timescales of the 
fueling mechanism(s) that are relevant to Seyfert activity. 
Looking at the cold ISM directly using 0.5-1" resolution 
CO data, the NUGA team have assessed the potential for 
gravitational torques fueling AGN and find that only 1/3 
of a sample of 25 galaxies are driving gas into the central 
100 pc (Haan et al. 2009, Garcfa-Burillo, S. & Combes 
|2012). However, a large fraction of the galaxies in 
this sample are low luminosity AGN with LINER or 
transition object classifications, rather than exclusively 
Seyferts. It is therefore possible that many of the galax- 
ies in the NUGA sample can be driven by more subtle 
fueling mechanisms on smaller scales that do not repre- 
sent the fueling requirements of higher luminosity local 
and higher redshift AGN. 

To constrain what mechanisms are responsible for 
fueling of local AGN, and potentially AGN at 
higher redshifts, we have carried out a matched 
Seyfert-quiescent galaxy study in which we mea- 



sure the 2D kinematics of molecular hydrogen that 
is ex pected to trace t he cooler ISM fueling AGN 
(e.g. iHicks et al.' '2009", 'Storchi -Bergmann et all 120101 
iRiffel fc Storchi-Bergmann 2011). In this first of a series 
of papers, we present a global analysis of the properties 
of Seyferts compared to quiescent galaxies and identify 
several differences in the nuclear (within a radius of 250 
pc) stellar and molecular gas properties. Future papers 
will present a more detailed analysis of the kinematics 
in the observed galaxies, including modeling to reliably 
demonstrate and quantify the detected ISM inflow. Our 
observed sample, data processing techniques, and con- 
firmation of nuclear activity in the selected galaxies is 
discussed in § 2. In § 3 the distribution and kinematics 
of both the nuclear stars and molecular gas are presented. 
Results based on our comparative analysis of the global 
properties in the Seyfert and quiescent galaxy subsam- 
ples are discussed in § 4, and in § 5 evidence for a lack 
of bulk H2 outflows and subsample differences in the de- 
tection of ionized gas emission are considered. This is 
followed by an interpretation of the observed subsample 
differences in § 6, and conclusions are presented in § 7. 

2. SAMPLE SELECTION, OBSERVATIONS, AND DATA 
PROCESSING 

2.1. Matched Seyfert- Quiescent Galaxy Sample 

In order to identify differences between Seyfert and 
quiescent galaxies that can provide insight into what 
mechanisms are fueling the Seyfert AGN the observed 
sample was selected to reduce any potential biases due 
to large-scale (>kpc) host galaxy properties. The sam- 
ple is therefo r e draw from the morphological study of 
IMartini et al.l (|2003al lH) who matched a sample of 26 
AGN and inactive galaxies pairs imaged with HST based 
on the host galaxy's Hubble type, i?-band luminosity, he- 
liocentric velocity, inclination, and angular size (see Ta- 
ble[T]for matching criteria). From this sample, pairs have 
been selected in which both galaxies are visible from VLT 
and have <200pc/" (distance < 40 Mpc), and for which 
the AG N is reported to be a Seyfert in IMartini et al.l 
(12003a). To ensure that differences in the sub-samples 
are linked to AGN activity level, weaker levels of nuclear 
AGN activity (i.e. tho se galaxies classified as LINERs 
or transition galaxies in IMartini et al.ll2003aD have been 
avoided to maintain a clear distinction between activity 
resulting from fueling of a BH and more global activity, 
such as vigorous star formation. The observed sample is 
comprised of 10 galaxies, five Seyfert and five quiescent 
galaxies, the host galaxy properties for which are given 
in Table [1] It is worth noting that such a matched sam- 
ple is not strictly necessary since, as discussed in 21 no 
large-scale properties have actually been found to corre- 
late with Seyfert AGN and therefore any differences are 
not expected to actually bias the samples. The observed 
sample is therefore also treated as two subsamples of five 
active and five quiescent galaxies, in addition to consid- 
ering a matched pair-wise comparison. 

For the specific 5 Seyfert-quiescent galaxy pairs ob- 
served there are no statistically significant systematic dif- 
ferences in any of the five parameters on which the pairs 
were matched (e.g. there is no systematic difference in 
inclination; see Table [Ij . //-band luminosities and disk 
sizes were also compared for the pairs, and the subsam- 
ples as a whole, and some marginal differences are found. 
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Table 1 

Galaxy Properties 



I.D. 


Galaxy 


AGN 


T 


Mb 


V 


R23 


/™C20 


Nuc. Dust 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


la 


NGC 3227 


Sey 1.5 


1.0 


-20.4 


1152 


0.17 


0.06 


C 


Iq 


IC 5267 


no 


0.0 


-20.3 


1713 


0.13 


0.06 


C 


2a 


NGC 5643 


Sey 2 


5.0 


-20.9 


1199 


0.06 


0.07 


GD 


2q 


NGC 4030 


no 


4.0 


-20.9 


1460 


0.14 


0.16 


TW 


3a 


NGC 6300 


Sey 2 


3.0 


-20.6 


1110 


0.18 


0.07 


C 


3q 


NGC 3368 


no 


2.0 


-20.3 


897 


0.16 


0.04 


cs 


4a 


NGC 6814 


Sey 1.5 


4.0 


-20.5 


1563 


0.03 


0.11 


GD 


4q 


NGC 628 


no 


5.0 


-20.2 


656 


0.04 


0.03 


N 


5a 


NGC 7743 


Sey 2 


-1.0 


-19.4 


1658 


0.07 


0.11 


LW 


5q 


NGC 357 


no 


0.0 


-19.9 


2541 


0.14 


0.14 


N 


Seytert Subsample 




2.4 


-20.4 


1336 


0.10 


0.05 










±2.4 


±0.6 


±254 


±0.07 


±0.02 




Quiescent Subsample 


2.2 


-20.3 


1453 


0.12 


0.10 










±2.3 


±0.4 


±741 


±0.05 


±0.04 





Note. — Galaxy I.D.s given in column 1, and used throughout the paper, 
are based on the matched pairs as described in ^ \2.1\ for example the first pair 
is the "active" Seyfert galaxy la and the "qui escent" galaxy Iq. Ifo st galaxy 
properties in columns 4-8 are as reported in IMartini et al.l ()2003al ) and the 
respective matching criteria for galaxy pairs are as follows: Hubble Type AT 
< 1, blue luminosity AMb < 0.5 mag, heliocentric velocity Av < 1000 km 
s~^, inclination Ai?25 < 0.1 (defined as the mean decimal logarithm of the 
ratio of the major isophotal diameter to the minor isophotal diameter), and 
the fraction of the angular diameter D25 qnhc host galaxy within the 20" field 
of view of the HST maps published in Martini ct al. 2003a A frac20 < 0.1. 
Column 9 gives the nuclear dust structure as classified by M artini e t al. 2003<i 
where the classification codes are GD for grand-design nuclear spiral, TW 
for tightly wound spiral, LW for loosely wound spiral, CS for chaotic spiral, 
C for chaotic circumnuclear dust, and N for no circumnuclear dust structure. 
Average values for each subsample are given in the bottom rows with the 
mean in the top row and the standard deviation in the row immediately 
below. 



The i7-band magnitudes, which can be taken to be rep- 
resentative of galaxy mass, are found to be about 24% 
(0.3 magnitudes; Table[2]) fainter on average in the active 
galaxies. If the pair containing NGC 628 is excluded, as 
is justified in !j3.2[ then this difference drops to Seyferts 
being only 13% fainter. In addition, the galaxy sizes are 
about 23% smaller in the active galaxy sample compared 
to the quiescent galaxies (Table [5]). In both cases the dif- 
ference is only Ict or less, and it could explain as much 
as, but not more than, an approximately 25% systematic 
difference in the subsamples due to the size and / or lumi- 
nosity based on the sample selection. Any difference in 
the individual bulge and disk components of the galaxy 
pairs is also likely to be limited to a similarly small differ- 
ence since the bulge-to-total flux ratio is closely related 
to Hubble type (iLaur ikainen et a"Lll2007f) , which was used 
as a criterion in matching the galaxy pairs (see Table 
[2]). A comparison of the large-scale bar classification, 
which was not considered in the pair matching, shows 
that bars are found in all five Seyferts and two of the qui- 
escent galaxies (NGC 3368 and NGC 357; IMartini et all 
l2003aD . There is therefore also no significant subsam- 
ple bias due to the presence of a large-scale bar, which 
is consistent with the lack of a correlation between such 
bars and Seyfert act ivity found in larger samples (e.g. 
iComeron et alll200l ). 

2.2. Observations and Reduction of Integral Field Data 



Near inf rared integral field dat a was obtained with 
SINFONI (lEisenhauer et al.l IIM IBonnet et all [2001 
on the ESQ Very Large Telescope (VLT). The pixel 
scale of the R~4300 if-band spectra (approximately 1.95- 
2.45/.tm) is 0". 125x0". 25 and typically the field of view 
(FOV) measured is 8"x8". All observations were ob- 
tained in seeing-limited mode with mean point spread 
function (PSF) full-width-half-maximum (FWHM) of 
0.55±0.05" or 54±24 pc. This excludes NGC 4030, for 
which all data is of worse quality due to rapidly degrad- 
ing seeing conditions and is estimated to be around l"(or 
132 pc) based on the ESO seeing monitor, and may par- 
tially explain the lack of II2 detection. Integration times 
range from 50-140 minutes per galaxy. See Table [3] for 
the details of each dataset. 

The spatial resolution is determined from PSF calibra- 
tion frames taken immediately before or after the sci- 
ence frames. For a few of the Seyfert galaxies the PSF 
can also be estimated from the data themselves. This is 
achieved by applying the fact that the intrinsic equivalent 
width of the C O 2-0 ban dhead is expected to fall within a 
limited range (lOliva et a l. 1995, Forstcr Schrcibcr 200dl 
iDavies et al.l 120071 ). The observed equivalent width can 
therefore be used to determine the dilution of the stel- 
lar feature due to the AGN continuum emission (see 
§2.31 and §3.11 for details), which, at this resolution, is a 
point source (predict e d sizes are < 2 pc;|ja ffc ct al. 200|i 
iTristram et all [20071 IRaban et al.l[2009r PSF FWHM 
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Table 2 

Galaxy Properties 



I.D. Galaxy AGN L2_iofceV //-band (iso.) H-hand (total) O.5D25 B/T 

(105 L0) (mag) (mag) (") (kpc) 

(1) (2) (3) (4) (5) (6) (7) (8) (9) 



la NGC 3227 


Sey 1.5 


4141 


8.10 


7.93 


161 


16.4 


0.28 


Iq IC 5267 


no 




7.85 


7.68 


157 


23.1 


0.32 


2a NGC 5643 


Sey 2 


116 


7.59 


7.46 


137 


11.2 


0.08 


2q NGC 4030 


no 




7.68 


7.60 


125 


16.5 


0.10 


3a NGC 6300 


Sey 2 


4272 


7.29 


7.20 


134 


11.1 


0.14 


3q NGC 3368 


no 


5 


6.65 


6.57 


228 


11.6 


0.25 


4a NGC 6814 


Sey 1.5 


4889 


8.10 


7.95 


91 


10.1 


0.10 


4q NGC 628 


no 


3 


7.42 


7.04 


314 


15.1 


0.08 


5a NGC 7743 


Sey 2 


1310 


8.84 


8.66 


91 


8.4 


0.26 


5q NGC 357 


no 




8.91 


8.82 


72 


11.2 


0.32 


Seyfert Subsample 




2946 


7.98/7.96 


7.84/7.81 


123/131 


11.5/11.8 


0.17/0.19 






±2100 


±0.59/0.68 


±0.56/0.64 


±31/29 


±3.0/3.3 


±0.09/0.10 


Quiescent Subsample 




7.70/7.77 


7.54/7.67 


179/146 


15.5/15.6 


0.21/0.25 








±0.82/0.93 


±0.84/0.92 


±94/65 


±4.8/5.6 


±0.12/0.10 


Mean pair-wise difference 




0.28/0.18 


0.30/0.15 


-57/-15 


-4.0/-3.8 


-0.04/-0.06 








±0.37/0.34 


±0.47/0.37 


±104/53 


±2.4/2.7 


±0.05/0.04 



Note. — Columns 1-3 give galax y I.D.s, names, and class ification. The X-r ay 2-10 kcV luminosities given in column 
4 ar c from the l i teratu re: fla) IVasudevan et al. f201dl) , (2 a) iLutz ct al l 1^2004), ("3a) Bcckmann ct al. ( .2006.) , {3q) and 
f4q) |Ueda et al.| | |2001D , f4a) |Turner & Pound j U989D , ("5a) IPanessa et al. (2006). The lumi nositi es given for la and 5 a 
are corrected for these sources being Compton thick as discussed in ^^s udcvan c t al.| l |201(]D and iPanessa et al.l |[200^) 
respectively. Columns 5 and 6 are //-band magnitudes as reported in the 2MASS Extended Source Catalog~ J^^mt^^^al1 
I200a) , where they are the magnitude within an elliptical aperture defined by the Ks-hand 20 mag arcsec ^ isophote and 
the total magnitude as derived from the isophotal aperture plus the integration of the surface brightness profile that 
extends from the isophotal aperture out to about 4 disk scale lengths, respectively. C olumns 7 and 8 are t he galaxy 
radius, which is half the diameter D25 of the B— 25 mag arcsec"^ isophote as reported in lMartini et all | |2003£j) , and the 
distances used to trans late to units of kpc are those given in Table [3] The bulgc-to-total fiux ratios (B/T) estimated 
based on Hubble type ijLaurikainen et al.l|2007l) . which arc given in Table are listed in column 9. Average values of 
each these quantities for the two subsample are given in the bottom rows with the mean, e.g. < 0.5Z)25, actiue >i in the 
top row and the standard deviation in the row immediately below. The means and standard deviations of the pair-wise 
differences, e.g. < 0. 5X^25, actiue-0. 5X^25, guiesceTit >j for each quantity are given in the last two rows. In both cases the 
first value is based on all 5 galaxy pairs and the second is excluding the fourth pair with NGC 628, which was excluded 
from the sample when considering stellar properties (see ^3.2t . 



Table 3 

Summary of Observations 



.D. 


Galaxy 


D 


Ref.'' 


pc/" 


Tint 


PSF FWHM 






(Mpc) 






(min) 


(") 


(pc) 


la 


NGC 3227 


21.1 


1 


102 


50 


0.55 


56 


iq 


IC 5267 


30.3 


2 


147 


140 


0.61 


90 


2a 


NGC 5643 


16.9 


3 


82 


140 


0.49 


40 


2q 


NGC 4030=^ 


27.2 


4 


132 


50" 


0.66=' 


87'' 


3a 


NGC 6300 


17.1 


4 


83 


140 


0.48 


40 


3q 


NGC 3368 


10.5 


5 


51 


40 


0.58 


30 


4a 


NGC 6814 


22.8 


3 


111 


140 


0.51 


57 


4q 


NGC 628 


9.9 


6 


48 


100 


0.59 


28 


5a 


NGC 7743 


19.2 


7 


93 


140 


0.54 


50 


5q 


NGC 357 


32.1 


3 


156 


90 


0.62 


97 



Observing conditions degraded rapidly during the acquisition 
of the data for NGC 4030 and this PSF FWHM measurement 
reflects a lower limit. The actual PSF FWHM during the major- 
ity of the NGC 4030 data is estimated to be closer to 1" based 

on ESQ seeing monitor at Paranal. 

fa References for dista nces: (1) ITerrv et all [2003 [2) 

iWillick et al. 1997, (3) [Tulhl [l988l. (4) [Tally et al.1 [2009I . 
(51 .Freedman et al.i |200I1 . (6) lOlivares E. eFall [20Tol . (7) 
iJensen et al.ll2003l . 

values determined directly from the data are consistent 
with those determined from the calibration frames to 
within 0.02", which is consistent with expected variation 
in the ambient seeing with time between the acquisition 
of the data and calibration frames. 



The data were reduced using the SINFONI custom re- 
duction package SPRED (jAbuter et al.ll2006D . which in- 
cludes all of the typical reduction steps applied to near 
infrared spectra with the additional routines necessary 
to reconstruct the data cube. After background sub- 
traction, which is performed with sky frames obtained 
interspersed with the on-source exposures, the data are 
flat-fielded and corrected for dead/hot pixels. Telluric 
correction and flux calibration, which is accurate to 10%, 
are carried out using A- and B-type stars. Residuals from 
the OH lin e emiss i on are minimized using the methods 
outlined in iDaviesI ()2007D . 

2.3. Confirmation of AGN classification 

The luminosity of the non-stellar AGN continuum is 
estimated in each galaxy by measuring the equivalent 
width (EQW) of the 12CO(2,0) A2.2935 bandhead. The 
intrinsic equivalent width of this absorption is constant 
for a broad range of star formation histories and ages 
and a deviation by more than 20% of the contant value 
of 12 A, i.e. lowe r than 9.6 A, implies dilution by 
the AGN continuum (lOliva et al.lll995llF5rster Schreibeii 
I200a IDavies et al.ll2007D . Theabsence of such dilution 
does not rule out the possibility of an AGN buried in 
a very optically thick material (i.e. a Compton thick 
AGN) , but the presence of dilution can be taken as verifi- 
cation of an AGN continuum. Alternatively a very young 
(<10 Myr) stellar population can cause similar dilution, 
but there is no indication of such a population in our 
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Figure 1. Azimuthal average of the CO equivalent widtii. Signif- 
icant dilution is seen in only three galaxies, all Seyferts. Triangles 
are Seyfert galaxies while circles are quiescent galaxies, and all 
curves are labeled according to those defined in Table [Tl The hori- 
zontal dashed line represents the lower bound on an undiluted CO 
bandhead equivalent width (see 12.31 for details). A typical er- 
ror bar, representative of the standard deviation of the equivalent 
width of all pixels at each annulus, is shown on the right. 

sample galaxies. 

As shown in Fig. [1] there is significant dilution of 
the CO equivalent width in three of the Seyfert galax- 
ies: NGC 3227 (I.D. la as defined in Table ^, NGC 
6300 (3a), and NGC 6814 (4a). A 2D map of the 
AGN continuum, calculated assuming the non-stellar 
continuum is equal to the total continuum times (1 - 
EQW observed /EQW intrinsic), is showu in Fig. [2] for cach 
of these galaxies. There is no indication of significant di- 
lution from an AGN continuum in the remaining two ac- 
tive galaxies, NGC 5643 (2a) and NGC 7743 (5a), despite 
their classif ication based on opt ical line ratios as Sey fert 
2 galaxies (|Phillips et all 119831 and lHo et al.|[l997al re- 
spectively). However, there is a weak downturn in the 
CO equivalent width of NGC 5643 (2a), which reaches 
the 9. 6 A cut off for AGN dilution at the very center. 
The lack of significant dilution of the CO bandheads in 
these two galaxies suggests that they have lower levels of 
AGN activity in comparison to the others and/or they 
are Compton thick. In both cases there is evidence of 
their being Compton thick (as discussed below), and for 
NGC 7743 some studies fi nd that it is consistent with a 
low luminosity AGN (e.g. iGonzalez-Martm et aLll2009[ ). 
As will be discussed in 21 these two galaxies still exhibit 
nuclear molecular gas and stellar properties that set them 
apart from their paired quiescent galaxies following the 
same trends seen in the other Seyfert galaxies. In all 
quiescent galaxies there is no indication of an AGN con- 
tinuum, confirming their non-AGN classification with a 
relatively constant CO equivalent width across the field 
of view, including at the center as defined by the peak of 
the Tf-band luminosity. As will be discussed in NGC 
628 has a CO equivalent width significantly lower than 
all other galaxies in the sample, which can be explained 
by a relatively small contribution of older stars. 

The AGN contribution in each galaxy is also assessed 
by gathering available 2-10 keV X-ray luminosities from 
the literature (Table [2]). With the exception of Compton 
thick sources, AGN are expected to produce a significant 
level of 2-10 keV emission and it can thus be used as an 
indicator of such activity. Based on available measure- 
ments, it is found that the selected active galaxies have 
2-3 orders of magnitude higher X-ray luminosities com- 
pared to the quiescent galaxies in the sample (three of the 
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Figure 2. Maps of non-stellar continuum in three galaxies (NGC 
3227, NGC 6300, and NGC 6814) where significant dilution of the 
CO bandheads is seen due to an AGN continuum. The central 
circle represents the radius at which the dilution factor is 2. In all 
maps North is up and East is to the left. 

quiescent galaxies were not detected in surveys in which 
they were included). This difference in X-ray luminosi- 
ties suggests that our selected AGN are indeed signifi- 
cantly more active than the selected quiescent galaxies, 
and are consistent with a Seyfert classification. There is 
evidence of significant obscuration in three of the Seyfert 
:alaxies: NGC 3 227 (Vasudcvan ct al. 20 lJ), NGC 564 3 
Lutz et al.ll2004l iGuainazzi c t al. 2004 , Shu et al.ll2007[) . 
and NGC 7743 ()Panessa et a l. 2006). For two of these 
Compton thick sources, NGC 3227 and NGC 7743, an 
estimate of the intrinsic 2-10 keV emission (corrected for 
obscuration) is available and listed in Table [21 Also of 
note is that the three galaxies with the strongest non- 
stellar AGN continua are also those with the highest X- 
ray luminosities, that NGC 5643 (2a), which has only 
a hint of AGN continuum based on dilution of the CO 
bandheads, has moderate observed X-ray luminosities 
and is known to be highly obscured, and that NGC 
7743 (5a), which shows no indication of an AGN con- 
tinuum, has an observed X-ray luminosity comparable 
to th e brightest quiescent galaxy ()Gonzalez-Martm et al.l 
l2009t i.e. before correction for obscuration). 

3. CHARACTERIZATION OF NUCLEAR MOLECULAR GAS 
AND STARS 

3.1. 2D Distribution and Kinematics 

In all 10 galaxies the first two CO bandheads at 
A2.2935, 2.3226 /xm are detected well enough to reli- 



: NGC 628 (4q) 
TNGC 6814 (4o) 
: NGC 7743 (5o) 
TNGC 357 (5q) 
; NGC 6300 (3a) 
TNGC 3368 (3q) 

NGC 5643 (2o) 
NGC 4030 (2q) 
NGC 3227 (la) 
IC 5267 (Iq) 
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Figure 3. Spectra integrated over a 8"x8" field of view arbitrarily 
shifted for clarity. The summed spaxels for each galaxy have been 
corrected to remove the rotational velocity based on a smoothed 
mapped of the stellar velocity field and each spectrum has been 
deredshifted. 
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Figure 4. Pair 1: NGC 3227 (la) 2D maps of CO absorption (stars) and H2 emission (molecular gas). Maps are as labeled from left to 
right: observed flux, velocity, and velocity dispersion. 2D maps of CO are in the top row and H2 in the bottom row, as indicated in the 
lower right of each set along with the galaxy name and pair identification. In all maps north is up and east is to the left. 
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Figure 5. Pair 1: IC 5267 (Iq) 2D maps. See Fig. Hfor details. 

ably measure the velocity and velocity dispersion across 
a large fraction of the 8"x8" field of view. H2 1- 
S(l) A2.1218 fj,m emission is detected in 7 of the galax- 
ies, with all 3 of those with non-detections being quies- 
cent galaxies. Two of these non-detections are galaxies 
for which Martini et al. (2003a) find no detectable nu- 
clear dust structure, which is known to correlate with 
low molecular gas content (e.g. I Young fc Scovillelll99lL 
iDumas et al.ll2010t ). Spectra of each galaxy integrated 
over the measured field-of-view are shown in Fig. [3] 
The 2D flux distribution, velocity, and velocity dis- 



persion (cr) of the stellar and gas kinematics have been 
extracted by fitting the CO bandheads and H2 emission, 
respectively. Before extraction of the kinematics the data 
cubes were spatially binned to a minimum SNR with the 
Voronoi binning technique of CappcUari & Copin (2003J . 
This was done to increase the SNR in outer radii where 
otherwise the signal would be too low for reliable line 
profile characterization. A target SNR of 10 and 5 per 
bin was used for all stellar and molecular gas maps, re- 
spectively. This level of binning maintained as high a 
spatial resolution as possible within the central 200 pc 
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Figure 6. Pair 2: NGC 5643 (2a) 2D maps. See Fig. Hfor details. 
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Figure 7. Pair 2: NGC 4030 {2q) 2D maps. See Fig. Hfor details. 

while still improving the SNR in outer radii. 

The line profile of the H2 e mission is extr a cted using 
our code LINEFIT described in lDavies et"all (|2011f ) (see 
also iForster Schreibeii 120091 ) . Briefly, this code fits the 
emission line in the spectrum of each spatial element, or 
spaxel, with an unresolved line profile (e.g. a sky line) 
convolved with a Gaussian, as well as a linear function 
to the continuum. Similarly, the A2.2935, 2.3226 /im CO 
bandheads were simultaneously fit with a broadened stel- 
lar template at each spaxel (except in NGC 628 where 
only the first CO bandhead was fit due to the limited 
spatial region over which the signal-to-noise ratio (SNR) 
of the second bandhead is adequate for a robust measure- 
ment). Two late- type stellar standards were obtained for 
this purpose as part of our observing program and both 
are found to be well matched with our galaxy spectra. 
For the Seyfert galaxies a region just blueward of the sec- 
ond CO bandhead (from 2.318-2.323 fim) was excluded 
from the fit to avoid any contamination from the [Ca 
VIII] A2.3213 fim coronal line which is present in at least 
4 of the 5 active galaxies (see § 15. 2p . All reported velocity 
dispersion values are intrinsic values, i.e. an instrumental 
broadening of 35 km has been accounted for. The 
uncertainty of fits to both the CO bandheads and the 



II2 emission are estimated using Monte Carlo techniques 
by refitting with added noise of the same statistics as the 
data. This is done 100 times, and the standard deviations 
of the best-fit parameters are used as the uncertainties. 
Typically the uncertainties in velocity and velocity dis- 
persion for both the stars and gas are on the order of 
5-15 km s^^. Bins with poor quality fits (generally SNR 
below 2 or kinematic errors greater than 30 km/s) were 
also masked. It is on these masked binned maps that 
all subsequent results are based for both the CO band- 
heads and II2. The resulting 2D maps of the stellar and 
molecular gas distribution and kinematics are shown in 
Figures [4] through [131 

The stellar luminosities are corrected for AGN con- 
tamination based on the dilution of the CO absorption 
equivalent width. As discussed in ? I2.31 to within 20% the 
intrinsic K-band CO equivalent width is 12 A, and a de- 
viation from this value represents di lution from th e non- 
stellar (e.g. AGN) continuum (e.g. iDavies "eFaII [2007). 
The stellar fraction of the continuum is therefore equiva- 
lent to the ratio of the observed CO equivalent width to 
the intrinsic value, and all stellar luminosities reported 
are corrected using this estimate. 

The kinematic center of the galaxy is taken to be that 
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Figure 8. Pair 3: NGC 6300 (Sa) 2D maps. See Fig. Hfor details. 
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Figure 9. Pair 3: NGC 3368 (3q) 2D maps. See Fig. |4]for details. Note that for tliis particular galaxy the spatial coverage for H2 is 
dictated by the SNR which was only high enough for reliable emission line characterization in the central overlapping region of the four 
science frames obtained. The SNR of the CO bandheads was high enough that line characterization out to the edge of each science frame 
was possible and therefore wider spatial coverage was achieved. 



determined from doing a 2D Gaussian fit to the non- 
stellar (i.e. AGN) continumxi, when available, or the K- 
band continuum, which is dominated in all galaxies by 
the stellar continuum. For the three galaxies with strong 
non-stellar continua, the centers based on the non-stellar 
continua and the total if-band continua are consistent 
within the measurement error of a tenth of a pixel. 

3.2. The Unique Case of NGC 628 

Compared to the rest of the observed sample the stel- 
lar attributes of NGC 628 are unique in a number of 



critical ways that warrants excluding this galaxy from 
our comparative analysis of the stellar properties (but 
not the H2 properties) in our two subsamples. This is 
justified by the fact that in this galaxy: (1) the spec- 
tral slope is red in contrast to the blue slope seen in all 
other galaxies (see Fig. [3]), (2) its stellar luminosity is 
more than 40 times lower than the other galaxies within 
the field of view measured (out to a radius of 150 pc), 
(3) the CO bandhead velocity dispersion is 30 km s~^ or 
less throughout the field of view measured, which is 4.2 
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Figure 10. Pair 4: NGC 6814 (4a) 2D maps. See Fig. |4]for details. 
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Figure 11. Pair 4: NGC 628 (4q) 2D maps. See Fig. Hfor details. 

times lower than the rest of the sample (5.0 and 3.5 times 
lower than the quiescent and Seyferts galaxy subsamples, 
respectively), and (4) the equivalent width of the CO 
bandheads is around 5A with a peak value at the center 
of 6.8A, which is well below that expected if the stellar 
light is dominated by the late- type stars typically present 
on the 0.1-1 kpc scales observed. Each of these measure- 
ments is unchanged when considering a fit to the first 
two CO bandheads rather than just the first bandhead 
on which our presented results are based as discussed in § 
13.11 The low stellar velocity dispersion, equivalent width, 
and luminosity could suggest that the contribution from 
the bulge, at least on the Q.l-lkpc scales observed, is very 
weak or absent such that the near-IR continuum is not 
dominated by the typical late-type stars of a "classical" 
bulge. In this scenario, the very low stellar luminosity 
can also explain the redness of the spectral slope since 
even a small amount of warm/hot dust emission, which 
could be (and likely is) present in any of the other sam- 
ple galaxies, will have a much more significant impact 
on the total continuum of NGC 628. Any dust emission 
would dilute the CO bandheads produced by the rela- 
tively few late-type stars present on the observed scales 
in this galaxy. In addition, to heat the dust on spatially 



extended (^100 pc) scales, one requires a population of 
young (<10 Myr) stars. The near-IR continuum of these 
young stars, while weak, would further reduce the depth 
of the CO bandheads 

To test the above hypothesis the integrated spectrum 
of NGC 628 is fit with a model composed of a compos- 
ite template stellar spectrum and a dust component as- 
suming a simple black body spectrum (i.e. emissivity 
Qv—1). As shown in Fig. [T31 a fit to the spectrum 
of NGC 628 is consistent with very little stellar contin- 
uum generated by late-type stars and is instead domi- 
nated by a contribution from warm/hot dust emission. 
The best-fit stellar template, which was generated from 
a stellar library (with no reddening applied), consisted 
of a combination of only late-type supergiants and no 
F- or G-type stars. The best-fit dust temperature was 
370 K, which results in a stellar fraction of 0.33 at 2.29 
/im. The observed low CO bandhead equivalent width 
of 5A can also be used to roughly estimate the stellar 
fraction since, as discussed, its value should intrinsically 
be about 12A. This approach gives a stellar fraction of 
0.42 (=5A/12A.), which is reasonably consistent with the 
fraction derived from the spectral fit. The spectrum of 
NGC 628 is thus dominated by warm/hot dust emission 
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Figure 12. Pair 5: NGC 7743 (5a) 2D maps. See Fig. |4]for details. 
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Figure 13. Pair 5: NGC 357 (5q) 2D maps. See Fig. Hfor details. 

(potentially heated by a young stellar population which 
does not contribute to the near-IR emission) because 
of the lack of a significant co ntribution from late-type 
stars. Supporting this picture, iZou et ahl (|2011[ ) find ev- 
idence for the previously presumed spheroid to actually 
be composed of a disk-like stellar component, a nucleus, 
and nuclear spiral plus ring structures with young stellar 
populations. Interestingly, this conclusion suggests that 
the observed unique near-IR photometric and kinematic 
characteristics might be a useful way to identify galaxies 
for which the spheroid is a mis-classified disk-like struc- 
ture. 

4. DIFFERENCES IN GLOBAL PROPERTIES WITHIN A 
RADIUS OF 250 PC 

An analysis of the global properties of the matched 
sample shows some significant differences between the 
subsamples in properties that are likely related to fuel- 
ing of BHs in Seyfert galaxies. All observed properties 
derived from the line profile analyses described in § [3] 
(fiux distribution, velocity, and velocity dispersion) were 
analyzed in both the stars and molecular gas to identify 
differences between the Seyfert and quiescent galaxy sub- 
samples. This was carried out by considering a range of 
apertures and annuli, as shown in Figures fT5ll20l Circular 
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Figure 14. Fit to the integrated spectrum of NGC 628. The 
fit (red curve) overlaid on the integrated spectrum (blue curve) 
consists of a sum of two components: late-type stars (solid black 
curve) and 370 K dust emission (dashed black curve; assumes an 
emissivity Q„=l). 

apertures were used since there is no systematic differ- 
ence in the subsample large-scale inclinations, and thus 
the effect of using circular versus elliptical apertures will 
be similar across the two subsamples. For the luminosi- 
ties reported the uncertainties are dominated by the 10% 
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Figure 15. Mean stellar luminosity as a function of radius (a) intrinsic values and (b) normalized to 250 pc. Error bars are the standard 
deviations of measurements within each aperture. Triangles represent Seyfert galaxies while circles are quiescent galaxies. The label on the 
curve for each galaxy is as given in Table [Tl 
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Figure 16. Total CO luminosity in apertures of radius (a) to 100 pc and (b) 100 to 250 pc and (c) the ratio of these two. Triangles 
represent Seyfert galaxies while circles are quiescent galaxies. Error bars represent the 10% systematic uncertainty of the flux calibration. 
The probability that the two subsamples are drawn from different parent populations based on the r=100-250 pc aperture is 93.8%. 



accuracy of the flux calibration. For the kinematics (e.g. 
velocity dispersion), mean values are an average of aU 
bins uniformly weighted (i.e. not luminosity weighted) 
and uncertainties are taken to be the standard deviation 
of values within the aperture considered. 

For those properties where a difference exists between 
the Seyfert and matched quiescent subsample, the prob- 
ability that the two subsamples are drawn from distinct 
populations is calculated based on a cumulative binomial 
distribution. This computes the probability that a given 
number of galaxies in one subsample will have a mea- 
sured quantity equal to or greater than the same quan- 
tity in the pairs of the other subsample. For example, if 
the H2 luminosity is found to be higher for Seyferts in at 
least 3 out of 4 galaxy pairs the probability that these 
two samples are drawn from different parent population 
is 68.8%, and if at least 4 out of 4 are higher in Seyferts 
the probability is 93.8%. 

Properties for which a significant difference exists be- 



tween the paired subsamples are discussed below. 
Stellar Luminosity — As shown in Figures [12] and [TCI 

beyond a radius of 100 pc all 4 Seyferts have lower lumi- 
nosities compared to their quiescent galaxy pairs (ex- 
cluding pair 4 with NGC 628 as discussed in § 13.21) 
while the luminosities are comparable within a radius 
of 100 pc. The average surface brightness (Fig. [15]) 
at a radius of 200 pc is found to be about two times 
lower in Seyferts compared to quiescent galaxies (a ra- 
tio of the two subsample means is < T^co, quiescent > 
I < Tico, active > = 1.9 whilc the mean of the pair- 
wise ratios is < T,cO,quiescent /'S'CO^active > = 2.1). 

The total stellar luminosity (Fig. [TO]) within 250 pc 
is about 1.5 times higher in quiescent galaxies (a ratio 
of the two subsample means is < Lco,quiescent > / < 
Lco,active > = 1-5 and the mean of the pair-wise ratios 
is < Lco,quiescent/Lco,active > = 1-5). The difference 
beyond 100 pc shown in Fig. [16b implies a 93.8% prob- 
ability that the two subsamples are drawn from different 
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Figure 17. Radially averaged kinematical properties of CO: (a) mean CO velocity dispersion as a function of radius and (b) inclination 
corrected rotational velocity. The CO velocity dispersion is the mean of pixels at each radius bin uniformly weight (rather t han luminosity 
weighted), and velocities have been corrected for inclination based on the large scale host galaxy inclination angles reported in'Mart ini et al.1 
^)03a). The negative velocities beyond 150 pc for IC 5267 (Iq), which are unphysical, are consistent within the error of the measurement 
with no detected rotation in this galaxy. The larger symbols and solid curve represent the mean of each subsample and the errors bars are 
the standard deviation of the mean at that radius and the radial bin size. Triangles represent Seyfert galaxies while circles are quiescent 
galaxies. The label on the curve for each galaxy is as given in Table [T] The typical errors of both the dispersion and velocity measurements 
are 5-15 km s~^. 
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Figure 18. Mean uniformly weighted CO velocity dispersion of all pixels within apertures of radius (a) to 100 pc and (b) 100 to 250 pc 
and (c) the ratio of these two. Triangles represent Seyfert galaxies while circles are quiescent galaxies. Error bars represent the standard 
deviation of the dispersion measured in the pixels within the aperture considered. The probability that the two subsamples are drawn from 
different parent populations based on both the r=0-100 pc and r=100-250 pc apertures is 93.8%. 



populations. However, the difference in stellar luminos- 
ity disappears within a radius of lOOpc, where both the 
total stellar luminosities and average surface brightness 
are comparable. Normalizing the surface brightness pro- 
files at a radius of 250 pc (Fig. [TSb : see also Fig. [T6b ) 
shows that the stellar luminosity in Seyferts is more cen- 
trally concentrated in all Seyferts compared to all qui- 
escent galaxies. Of note is that the stellar luminosities 
in the Seyfert subsample have been corrected for AGN 
contamination based on the dilution of the CO bandhead 
assuming an intrinsic CO equivalent width of 12A (see § 
12.31 and § 13. II for details). Taking a more extreme intrin- 
sic CO equivalent width of 17A, whi c h is th e highest CO 
equivalent width in the lOliva et al.l (jl995f ) sample for a 
star-forming H II galaxy and lowers the central stellar lu- 
minosities by 30%, the Seyfert stellar luminosity profiles 
are still found to be more centrally concentrated than 
those of the quiescent galaxies. 

Stellar Dispersion — The stellar velocity dispersion 
also differs in Seyferts compared to quiescent galaxies, 
with all 4 pairs included in the analysis having a lower 
velocity dispersion out to a radius of 200 pc compared 
to their quiescent counterparts (see Figures flTk and fTS]) . 
Given these statistics, the probability that the Seyfert 



galaxies are drawn from a different parent population in 
which the velocity dispersion is lower compared to the 
quiescent galaxy population is thus 93.8%. The mean 
stellar dispersion in Seyferts out to a radius of 100 pc 
is about 1.4 times lower than that in quiescent galax- 
ies (< (TCO, quiescent > / < (JCO,active > — 1-4 and 

< CF c o. , quiescent /(y CO, active > = 1-5). Thcrc is an uptum 
in the velocity dispersion in all but one of the Seyfert 
galaxies at a radius of about 200-250 pc, which is not 
seen in the quiescent galaxy population, suggesting that 
the dispersion in Seyferts is comparable to that in quies- 
cent galaxies at radii shortly beyond this. 
H2 Luminosity — On the scales measured, the H2 lumi- 
nosity is observed to be elevated in all 5 Seyfert galaxies 
compared to their paired quiescent galaxies out to a ra- 
dius of at least 250 pc f Figures \W\ and 1^0)) . This compar- 
ison includes the 3-f7 upper limits derived for the three 
galaxies for which no H2 detection was observed. Fur- 
thermore, as shown in the normalized surface brightness 
profiles (Fig. [T9b). Seyferts have a steeper luminosity 
profile compared to the quiescent galaxies. On average 
Seyferts have 4 times higher H2 luminosity within a ra- 
dius of 100 pc compared to their quiescent counterparts, 
with a range in the sample up to 9 times higher luminos- 
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Figure 19. Mean H2 luminosity as a function of radius (a) intrinsic values and (b) normalized to 250 pc. Error bars are the standard 
deviations of measurements writhin each aperture. Triangles represent Seyfert galaxies while circles are quiescent galaxies. Dashed curves 
and arrows represent upper limits. The label on the curve for each galaxy is as given in Table [l] 
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Figure 20. Total H2 luminosity in an apertures of radius (a) to 100 pc and (b) 100 to 250 pc and (c) the ratio of these two. Triangles 
represent Seyfert galaxies while circles and upper limits are quiescent galaxies. Error bars represent the iO% systematic uncertainty of the 
flux calibration. The probability that the two subsamples are drawn from different parent populations based on the r=0-100 pc aperture 
is 96.9%. 



ity. The probability that the five Seyferts are drawn from 
a different parent population in which the H2 luminosity 
is higher compared to the population of quiescent galax- 
ies is 96.9% (and 93.8% at radii beyond 150pc where the 
limited field of view measured for NGC 628 limits the 
sample to 4 pairs). 

5. ASSESSMENT OF MOLECULAR GAS OUTFLOWS AND 
THE PRESENCE OF IONIZED GAS 

5.1. No evidence for Bulk Outflows in Molecular 
Hydrogen 

To determine if there is evidence of bulk outflows in 
the molecular hydrogen an independent fitting of the red- 
shifts of the H2 emission line and the stellar continuum 
have been performed with the integrated spectra shown 
in Fig. [3l Any systematic offset in the derived line and 
continuum rcdshifts would imply radial motion of the 
molecular gas with respect to the nuclear stellar popu- 
lation. This would reveal any outflows in the molecular 
gas, such as those often seen in ionized gas on the order 
of several hundreds of km s^^ or more, including in at 



least some of the Seyferts in this sample (e.g. NGC 6814; 
iMiiller-Sanchez et al.l [20TT). 

The H2 emission was fit in the integrated spectrum of 
each of the 7 galaxies in which it was detected, using a 
method analogous to that described in § 13.11 for fitting 
of the emission line at each spaxel in the 2D data cubes. 
The fit to the stellar continuum was similar to that de- 
scribed for NGC 628 in § 13.21 the rest frame wavelength 
region of 2.285-2.305 /im was fit with a stellar template 
generated from a stellar library (without a correction for 
reddening), but in this case assuming a fixed dust frac- 
tion and stellar dispersion. Although this fit is to the 
entire stellar continuum within the fit wavelength range, 
the resulting best-fit redshift is dominated by the CO 
bandheads. A comparison of the redshifts determined 
via each of these methods shows that there is no sys- 
tematic offset in velocity between the molecular gas and 
stellar continuum, and thus no detected bulk outflow of 
molecular gas. Any differences are at the level of about 5 
km s~^ (except NGC 5643 which is only 18 km s^^), and 
there is no obvious systematic trend in the Seyfert ver- 
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sus quiescent subsamples. This confirnis the conchisions 
of others that the kinematics of the nuclear molecular 
gas in Seyfert galaxies is dominated by rotation rather 
than bulk out flow and is thus a reasonable tracer of the 
cool ISM (e.g. [H icks et al. 2009, Stor chi-Bergmann et al.l 
[20T0I iRiffel fc Storchi-Bergniann.20ll 

5.2. General Properties of ionized Gas: Br 7 and 
Goronal Line Emission 

Brief descriptions of the ionized gas detected in the 
2D data cubes and integrated spectra shown in Fig. [3] 
are given below, including a comparison of the general 
properties of this emission in the Seyfert versus quies- 
cent subsamples. A detailed analysis of the nature of 
this emission, and its 2D distribution and kinematics, is 
beyond the scope of this paper, but will be presented in 
a future paper. 

Br 7 2.166 /im — The presence of Br 7 emission in a 
Seyfert galaxy can be interpreted as the result of star 
formation and/or the AGN activity itself (emission from 
the narrow and/or broad line regions). In the absence 
of other typical signatures of AGN activity (e.g. non- 
stellar continuum. X-ray emission), the detection of Br 
7 emission is expected to be the result of star formation. 
In the observed sample, none of the quiescent galaxies 
have measurable Br 7 emission in their integrated spec- 
tra or at any location in the 2D data cubes. This suggests 
that these galaxies have a relatively old stellar population 
within the nuclear region. On the other hand, all five of 
the Seyfert galaxies have detectable Br 7 emission. Fur- 
thermore, in each of these cases the 2D distribution of the 
emission extends beyond the FWHM of the PSF (with 
a SNR of at least 2, and much higher in several cases), 
and thus beyond the region significantly influenced by 
emission from the broad line region. The emission in 
the Seyfert galaxies could therefore be due to excitation 
within the narrow-line region, and thus associated with 
the AGN activity itself, but a contribution from rela- 
tively recent star formation in the nuclear region is also 
possible. A future paper will present 2D maps of the Br 
7 emission distribution and kinematics as well as explore 
the nature of this line emission. 

[Ca VIII] A 2.321 /im — Also of note is the presence 
of the AGN coronal line Ca [VIII] emission in at least 
4 of the 5 Seyfert galaxies, the exception being Seyfert 
NGC 7743. Any detection of this emission in NGC 7743 
is marginal, although this is not surprising since, as al- 
ready noted, there is evidence t hat t his is a weaker AGN 
than the others in the sample (? I2.3|) . In some cases, this 
emission is not obvious in the integrated spectra in Fig. 
[3] since it is generally concentrated to the vicinity of the 
AGN and is thus washed out in the full FOV integrated 
spectra shown. In none of the quiescent galaxies is this 
coronal line emission detected, providing additional con- 
firmation of the lack of AGN activity in these galaxies. 

6. DISCUSSION 

We have identified statistically significant differences 
between Seyfert and quiescent galaxies, and it is these 
properties unique to Seyfert galaxies that provide the 
first direct clues to what processes lead to fueling of 
their central BHs. We find that Seyferts have a steeper 
stellar surface brightness profile compared to quiescent 
galaxies and lower velocity dispersion out to a radius of 



200 pc. In addition, Seyferts have an elevated H2 lumi- 
nosity within this same nuclear region that is more cen- 
trally concentrated than the gas distribution seen in qui- 
escent galaxies. We therefore conclude that we are now 
pro bing scales relevant to the fueling of Seyfert galax- 
ies. iDumas et al. (201)7) also conclude that it is on scales 
of hundreds to tens of parsecs that a link between host 
galaxy and BH exists based on their analysis of disk of 
ionized gas which are likely the counterparts to the cold 
gas disks detected in H2. This suggests that the mecha- 
nisms responsible for the inflow of material are likely to 
be shaping the observed characteristics of the nuclear (r 
< 250 pc) stars and gas. By understanding how these 
observed differences arise, combined with future detailed 
modeling of the observed kinematics, we can begin to 
form a more complete picture of what processes are dic- 
tating the growth of BHs in the local universe and how 
they relate to coevolution of BH and host galaxy. 

6.1. Significance of Differences in Molecular Gas 
Properties 

If the H2 lumi nosity reflects gas mass, as it does on 
_large scale s (e.g. Dale et all 120051 IMiiller Sanchez et al.l 
20090 . then the active galaxies have gas 
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masses on average 4 times higher than similar quiescent 
galaxies within a radius of 100 pc. On the other hand, it 
is known that the H2 luminosity on small scales can be 
strongly influenced by local excitation. For example, the 
presence of shocks or a young stellar population can influ- 
ence the H2 emission. X-ray heating of the gas from the 
AGN itself is al so a possibili t y wh ic h has support both 
theoreticallv (iMaloney et al.l Il996l iMeiierink &: Spaajisl 
120051 iMeiier ink et al.l I2007D as well as obse rvationallv 
(|Malonevl ll997. Rodrig uez- Ardila et al.l[2005f ). Thus the 
enhanced H2 luminosity could reflect different excitation 
mechanisms and efficiency (i.e. fraction of hot gas) be- 
tween the two subsamples. In either case, the difference 
in H2 luminosity points to a physical difference between 
active and quiescent galaxies on scales of hundreds of 
parsecs. 

6.2. Significance of Differences in Stellar Properties 

The lower stellar velocity dispersion in Seyfert galaxies 
can be understood to be either an indication of a differ- 
ence in the orbital structure of the nuclear region (e.g. 
the stars are in a dynamically colder disk-like structure) 
or a difference in mass distribution (e.g. the nuclear re- 
gion measured contains less mass). The data presented 
here indicate that Seyfert and quiescent galaxies differ in 
the orbital structure of the nuclear region. A key piece 
of evidence leading to this conclusion is that the lower 
stellar velocity dispersion observed in Seyferts is not ac- 
companied by a lower rotational velocity. As can be seen 
in Fig. [T7b . at all radii measured there is no statistically 
significant difference in the rotational velocities of the 
paired galaxies or in a comparison of the subsamples as 
a whole. Thus a given rotational velocity is accompanied 
by a lower velocity dispersion in Seyferts than it is in qui- 
escent galaxies, which implies that the contribution from 
the dynamically hot stellar component is less significant 
in Seyferts than it is in quiescent galaxies. 

In an effort to establish the relative contribution of 
rotational and random velocities in the nuclear region, 
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Figure 21. Global effective CO V/cr, following |Binnevll2005l . (a) as a function of radius and (b) versus moment cUipticity e within a radius 
of 100 pc. Velocities has been corrected for inclination based on the large scale host galaxy inclination angles reported in lMartini et all 
l[2003al ) . Triangles represent Seyfert galaxies while circles are quiescent galaxies. The label on the curve for each galaxy is as given in Table 
11^ 
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Figure 22. Mean CO y/v'^ + a'^ as a function of radius. The CO velocity dispersion is the mean of pixels at each radius bin uniformly 
weight (rather tha n luminosity weighted ), and the velocity has been corrected for inclination based on the large scale host galaxy inclination 
angles reported in lMartini et al.1 l l2003al V The larger symbols and solid curve represent the mean of each subsample and the errors bars are 
the standard deviation of the mean at that radius and the radial bin size. Triangles represent Seyfert galaxies while circles are quiescent 
galaxies. The label on the curve for each galaxy is as given in Table[T] The typical errors of the \/v^~+~a^ quantity are 10-20 km s^^. 



we consider the effective V/a (a measure of dynami- 
cal state of t he system in terms of ordered versus ran- 
dom motions: lBinnevll2005[) . the specific angular momen- 
tum \fi (a measure of dynamical state in terms of the 
observed projected a ngular momentum per units mass; 
lEmsellem et al.lf2007| ). and the moment ellipticity e, all 
of which are generally independent of mass. Unlike the 
quantities discussed above, e.g. tr, which are local val- 
ues generally within an aperture, these are global val- 
ues within the stated radius. All velocities have been 
corrected for inclination using t he large-scale values re- 
ported in iMartini et all (|2003al ) (R25 in Table HJ. As 
the radial profiles of the effective V/a shown in Fig. [2Tk 
illustrate, at a radius of 100 pc, each of the four pairs 
has a higher value for the Seyfert galaxy (the pair in- 
cluding NGC 628 is excluded as justified in §[33). This 
is also shown in Fig. [2Tb where the effective V/a is 



plotted against e for the global values within a radius of 
100 pc. The result is the same for the specific angular 
momentum A^, which is expected given this quantity is 
related to V/a for regular velocity fields. At a radius of 
200 pc (Fig. [2Tk). three of the pairs still have higher 
values in the Seyferts, and Pair 3 has comparable val- 
ues (the Seyfert galaxy 3a is just a bit lower than the 
quiescent galaxy 3q). Beyond this radius the remaining 
three pairs for which a comparison can be made have 
higher values in the Seyferts. These quantities therefore 
indicate that, at least within a radius of 200 pc, the sub- 
samples differ in orbital structure, with Seyferts having 
a dynamically cooler structure compared to the quies- 
cent galaxies. Outside of a radius of 200 pc the trend for 
increasing stellar velocity dispersion with radius in the 
Seyfert galaxies suggests that the subsample differences 
in orbital structure disappear beyond this radius. 
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A lower velocity dispersion could also be the result of 
a systematic subsample difference in inclination, but we 
consider this unlikely due to (1) a lack of systematic dif- 
ference in inclination on large scales (see Table [1] mean of 
the difference in pairs is 7.5±5.3° lower in Seyferts) and 
(2) a lack of systematic difference in small-scale inclina- 
tion derived from fitting of the stellar continua isophots 
shown in Figures l4l ll3l (mean of the difference in pairs is 
7.8±5.5° higher in Seyferts). A lack of systematic dif- 
ference in the moment ellipticities (Fig. 121b ) provides 
additional support for the conclusion that it is not an 
issue of inclination that results in the observed differ- 
ences in velocity dispersion. In addition, the subsample 
differences in effective V j a and specific angular momen- 
tum Aij cannot be explained via a systematic difference 
in inclination since there is no systematic difference in 
the observed moment ellipticity e. It is therefore more 
likely that a difference in orbital structure explains the 
observed differences in V ja and rather than a sys- 
tematic difference in inclination. 

In order to assess the mass in the nuclear region, the 
quantity \Jv'^ + a"^ is considered (Fig. [22)) . This quan- 
tity gathers the first velocity moments so that a de- 
crease in the dispersion should be partly compensated 
by an increase in velocity if the difference in disper- 
sion were fully due to a y ariation in the enclosed mass 
(jBinnev fc Tremainell2008l) . As before, all velocities have 
been corrected for inclination using t he large-scale val- 
ues reported in iMartini et all ()2003aD (R25 in Table [1]). 
Throughout the region measured, the difference between 
the Seyfert-quiescent pairs is not statistically significant 
when considering y/v^ + a^. The only pair with a signif- 
icant difference in this quantity is Pair 5 in which the 
quiescent galaxy has a higher value throughout the ma- 
jority of the FOV measured. In the remaining pairs the 
difference between the Seyfert and quiescent galaxies is 
within the error of the measurement (a difference of less 
than 25 km s^^). There is also no significant difference 
in the subsample means (solid curves in Fig. [22)) . This 

lack of subsample difference in ^/^P + cJ^ supports the 
theory that the lower stellar velocity dispersion observed 
in Seyferts galaxies is not driven by a difference in mass 
enclosed within the nuclear region. 

In addition to kinematical considerations, a number 
of photometric quantities suggest that the lower stellar 
velocity dispersion in the Seyfert subsample can be at 
most only partially explained by a systematic difference 
in the subsample mass within the nuclear region. If the 
observed subsample difference in dispersion of a factor 
of 1.4 within a radius of 100 pc is due to a difference 
in mass, and if it is assumed that mass is a function of 
(7^, then a mass difference of a factor of two is required. 
The difference implied by the lower total stellar lumi- 
nosity in Seyferts within a radius of 250 pc (Fig. (THl see 
also discussion in Q , if taken to reflect mass and assum- 
ing a fixed mass-to-light ratio, is a mass difference of a 
factor of only 1.5±0.5. Considering the fl-band 2MASS 
estimated total photometry (Table [2]) as an indicator of 
stellar mass, the mean difference in the subsamples is 
0.15 ± 0.37 mag, with Seyferts slightly fainter (and thus 
less massive). This suggests a mass difference of a fac- 
tor of l.l5to.tt 

also less than the full difference of 2 
required to explain the lower stellar dispersion. A com- 



parison of the subsample galaxy sizes, quantified in this 
case as half of D25 (Table [2]) , reveals a subsample dif- 
ference with Seyferts being 23±13% smaller, and thus 
less massive, than the quiescent galaxies. However, a 
systematic difference in D25 or //-band photometry does 
not necessarily reflect a difference in actual nuclear mass 
since these quantities are signiflcantly influenced, if not 
dominated, by properties of the stellar disk. A compar- 
ison of the bulge-to-total flux ratios estimated based on 
the Hubble type (see Table [2] for details) indicates there 
is no statistically signiflcant difference between the sub- 
samples. Therefore the difference in disk mass implied 
by subsample differences in D25 and //-band photometry, 
as well as the difference indicated by the total stellar lu- 
minosities, suggest that the total nuclear masses may dif- 
fer, but not signiflcantly enough to explain the observed 
difference in velocity dispersion. In addition, the steeper 
gradient of the central stellar surface brightness observed 
in the Seyfert galaxies compared to the quiescent galaxies 
cannot be easily explained through a difference in mass 
and suggests a more complex situation. 

Having ruled out a systematic difference in inclina- 
tion and flnding minimal evidence for a subsample dif- 
ference in mass within the nuclear region, we conclude 
that the observed subsample difference in stellar velocity 
dispersion is due primarily to a difference in orbital struc- 
ture. In this scenario Seyfert galaxies have a dynamically 
cooler system within a radius of 200 pc in comparison to 
quiescent galaxies which results in the lower observed 
nuclear stellar velocity dispersions. 

6.3. Evidence for ISM Inflow 

The ways in which Seyfert galaxies are observed to dif- 
fer from quiescent galaxies, which presumably also har- 
bor central BHs that are in a dormant state, suggests 
that the fueling of an AGN is accompanied (or triggered) 
by an altered nuclear region that extends out to a radius 
of 200-250 pc. The observed lower stellar velocity dis- 
persion within this region appears to be due to a com- 
paratively dynamically cooler disk-like orbital structure 
in Seyfert galaxies. These so called "cr-drops", can be 
understood as evidence for the presence of an additional 
nuclear stellar component in Seyfert galaxies that is dy- 
namically cold in comparison to the (dynamically hot- 
ter) stars that dom inate the detected stellar light in the 
quiescent galaxies (Emsellem et al. ' 2001llWozniak et al.l 
,2003,Wozn iak & Cha m pavert 2006). Similar "a-drops" 
have been obse r yed in a number of AGN galaxies (e.g . 
Emselleml^OOSl IComeron et all [20081 IRiffel et al.l[2009l 
201 a 1201 llT and there is also evidence that a central de- 
crease in stell ar dispersio n is correlated with Seyfert nu- 
clear activity (jPumas et al. 2007, Comeron et al. 2008). 
This interpretation of the lower stellar velocity dispersion 
as a dynamically cold nuclear component is consistent 
with the more centrally concentrated stellar distribution 
seen in the Seyfert stellar light proflles: as the dynam- 
ically cold nuclear component begins to dominate the 
light within the central 200 pc the proflle represents that 
of the nuclear component rather than the higher velocity 
dispersion stars and thus steepens. In quiescent galaxies 
such a nuclear stellar structure may either be absent or 
be much harder to detect because its stellar population 
no longer dominates the near-IR light over the higher ve- 
locity dispersion component (e.g. the dynamically cold 
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stars are older and fainter in comparison to that found 
in Seyferts). 

A consequence of interpreting the observed differences 
in subsample stellar properties as evidence for a nuclear 
dynamically cold component in Seyferts is that a reser- 
voir of gas must be present from which the stellar pop- 
ulation of this component can form, and thus inflow of 
ISM must occur to supply the gas to the nuclear region. 
The presence of such a nuclear component thus requires 
that ISM inflow has occurred. In this scenario the ob- 
served enhanced H2 emission and greater central concen- 
tration in Seyferts can then naturally be interpreted as 
evidence of an enhanced nuclear gas mass rather than 
being explained solely by a difference in excitation or 
efficiency. The central BH can then be fueling either 
directly from this supply of nuclear gas or indirectly 
throug h smaller scale processes, including from s tellar 
winds (jDavies et al.l [20071 iSchartmann et all 120101 In 
such a scenario, the Br 7 emission, which was detected 
in only the Seyfert subsample, may be partially due to 
excitation from young/hot stars within this nuclear re- 
gion of enhanced gas mass and nuclear star formation. 

Also of note is that in galaxies from both subsaniples 
there is evidence of non-circular velocities (e.g. radial 
flow) in the velocity fields of II2, despite rotation gen- 
erally dominating the gas motion. This can be easily 
seen in Fig. l4llT3l through a comparison of the stellar 
zero iso- velocity contours with that of II2 , which is often 
twisted in the nuclear region as expected in the presence 
of a nuclear bar or spiral. Detailed modeling and char- 
acterization of the nuclear kinematics of each galaxy will 
be presented in a future paper and any additional sub- 
sample differences will be identified to determine what 
fueling mechanisms contribute to the fueling of the AGN 
present in the Seyfert subsample. 

7. CONCLUSIONS 

With a sample of matched Seyfert-quiescent galaxy 
pairs we have started to identify the processes that drive 
BH growth in the local universe, as well as the processes 
potentially at work at higher redshifts, by directly mea- 
suring the morphology and kinematics of the stars as well 
as the warm molecular gas that traces the ISM making 
up the reservoir of AGN fuel. We find subsample dif- 
ferences within the central kpc that are correlated with 
AGN activity. The subsample differences indicate that 
Seyferts have (1) a more centrally concentrated nuclear 
stellar surface brightness with a lower stellar luminos- 
ity outside of a radius of 100 pc, (2) a lower stellar ve- 
locity dispersion within a radius of 200 pc, (3) elevated 
H2 luminosity out to a radius of at least 250 pc, and (4) 
more centrally concentrated H2 surface brightness pro- 
files. Although this sample is rather small, comprising 
of 5 Seyfert-quiescent galaxy pairs, the subsample differ- 
ences are statistically significant. This suggests that with 
these observations of the central 1 kpc down to 50 pc we 
are now probing scales relevant to the link between the 
host galaxy and BH growth. In addition, Br 7 emission 
is detected in all 5 Seyfert galaxies and the coronal line 
[Ca VIII] is present in 4 of the 5 Seyferts. In contrast, 
both emission lines are undetectable in the spectra of 
all of the quiescent galaxies, confirming the lack of AGN 
activity in these galaxies. 

A likely interpretation of the observed subsample dif- 



ferences is that in Seyfert galaxies a nuclear component 
composed of both stars and gas is present that is either 
undetected or not present in quiescent galaxies. This is 
possible if the nuclear component dominates the detected 
stellar light out to a radius of 100-250 pc and is dynam- 
ically cold in comparison to stars detected in quiescent 
galaxies. The elevated H2 luminosities and higher cen- 
tral concentration would then represent a reservoir of gas 
mass available for forming the nuclear stellar component 
and potentially fueling the BH. Assuming the enhanced 
H2 is a refiection of gas mass (as opposed to being en- 
tirely due to a difference in excitation mechanisms), this 
then implies a factor of 4 difference on average in gas 
mass between Seyfert and quiescent galaxies within a ra- 
dius of 100 pc. Under this scenario the nuclear disk-like 
structure observed is generated by inflow to the nuclear 
region and likely fueling the BH either directly or indi- 
rectly through smaller-scale processes occurring within 
this nuclear disk. Another possible explanation of the 
observed subsample differences is that Seyferts have less 
mass within a radius of 250 pc compared to quiescent 
galaxies, however we find minimal evidence for this sce- 
nario. 

Future work will include a detailed analysis of the kine- 
matics observed in the sample and a subsample compari- 
son will further identify any additional properties unique 
to Seyfert galaxies. Modeling of the observed kinematics 
will then classify (e.g. identify nuclear spirals, bars) and 
quantify any on-going inflow of gas within the central 
kpc. With an understanding of how Seyferts differ in 
comparison to quiescent galaxies in their nuclear stellar 
and gas characteristics, and most critically in their inflow 
of gas, on scales from 500 down to 50 pc we can begin 
to identify what processes are essential to BH growth in 
local galaxies and then explore how this might extend to 
higher rcdshift systems and impact the general evolution 
of the host galaxy. 
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